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Аннотация 

Введение. При конечно-элементном моделировании высотных зданий из монолитного железобетона, как пра- 
вило, используется рамно-стержневая расчетная схема. Данная модель не позволяет получить данные о распре- 
делении компонент тензора напряжений в зоне сопряжения перекрытия с колонной. Поэтому проблема совер- 
шенствования методики прочностного расчета в местах стыков перекрытий и колонн является актуальной. Це- 
лью настоящего исследования является разработка концепции фрагментации рамного каркаса для оценки не- 
сущей способности перекрытий. На тестовом примере шестипролетного трехэтажного монолитного железобе- 
тонного каркаса выполнены численные эксперименты с использованием объемно-стержневой и объемно- 
пластинчатой моделей повторяющегося конструктивного фрагмента. Разработаны практические рекомендации 
для уточненного прочностного расчета перекрытий монолитных железобетонных каркасов многоэтажных 
зданий. 

Материалы и методы. Вычислительные эксперименты выполнены с помощью программного комплекса 
АМЗУ$ Меспашса], в котором реализован метод конечных элементов в форме метода перемещений. Для моде- 
лирования напряженно-деформированного состояния монолитного железобетонного каркаса применен пла- 
стинчато-стержневой ансамбль конечных элементов. Уточненный расчет зоны сопряжения перекрытия и ко- 
лонны при статическом нагружении выполнен с использованием ансамблирования объемных, балочных, фер- 
менных и пластинчатых элементов. 

Результаты исследования. Разработана инженерная методика численного анализа напряженно- 
деформированного состояния зоны сопряжения перекрытия и колонны железобетонного монолитного каркаса 
при статическом нагружении. Наиболее точный результат обеспечила конечно-элементная модель, построенная 
с использованием балочных конечных элементов в качестве арматурных стержней. 

Обсуждение и заключения. Разработанная методика численного моделирования сопряжения перекрытия и 
колонны позволила оценить реальный запас прочности данного узла с учетом реальной геометрии армирующих 
сеток, а также уточнить несущую способность монолитного железобетонного каркаса при различных сценариях 
нагружения. 


Ключевые слова: метод конечных элементов; объемные, балочные, ферменные, пластинчатые конечные эле- 
менты; безригельное перекрытие с бескапительным стыком; модель сопряжения плиты перекрытия и колонны; 
модели дискретного армирования перекрытия. 
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Введение. В настоящее время проектирование многоэтажных высотных зданий из монолитного железобе- 
тона базируется на рамно-связевой конструкционной схеме, позволяющей в известной степени обеспечить 
«живучесть» здания в случае прогрессирующего (лавинообразного) разрушения [1, 2]. Отметим, что в отече- 
ственной строительной науке под термином «прогрессирующее разрушение» понимается процесс обрушения 
несущих конструкций на нескольких этажах здания или на одном этаже площадью более 80 м”. Данное явление 
возникает в результате одномоментного разрушения, как правило, одного несущего элемента, с последующим 
быстро нарастающим обрушением всего здания или его части по сценарию «домино». 

В качестве основных конструкционных элементов монолитных железобетонных каркасов многоэтажных 
зданий выступают повторяющиеся фрагменты колонн и безригельных перекрытий, связанных бескапительны- 
ми стыками [1]. 

Несмотря на многолетний опыт проектирования зданий с рамно-связевыми каркасами, в мировой практике 
известны случаи прогрессирующего разрушения данных объектов. Причины таких явлений главным образом 
обусловлены ошибками при проектировании зоны сопряжения перекрытий и колонн в сочетании с нарушением 
установленных правил эксплуатации зданий. Целью настоящего исследования является разработка инженерной 
методики прочностного расчета рамного каркаса многоэтажного здания из монолитного железобетона с учетом 
объемного характера напряженного состояния в зоне стыка перекрытия и колонны. 

При конечно-элементном моделировании изгибаемых железобетонных конструкций обычно используют 
подход, основанный на представлении бетона двумерными или трехмерными конечными элементами (КЭ), по- 
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строение которых базируется на принципах теории упругости. Армирующие стержни, как правило, моделиру- 
ются балочными или ферменными КЭ соответствующей размерности. По способу ансамблирования объемных 
и стержневых КЭ различают следующие схемы представления армирования [2, 3]: дискретно-распределенное, 
при котором координаты узлов разнотипных элементов совпадают (рис. 1 а); армирование с использованием 
так называемых встроенных конечных элементов (рис. 1 6). В последнем случае координаты узлов разнотип- 
ных КЭ не совпадают и предусматривается процедура «конденсации» элементов матриц жесткости стержнево- 
го элемента на смежные узлы объемного (основного) элемента (рис. 1 в). Отметим, что вычислительная техно- 
логия встроенных КЭ применима исключительно для решения задачи об обобщенном плоском напряженном 
состоянии. Причем основные КЭ должны быть изопараметрическими и поликвадратичными, а встроенные 
КЭ — прямолинейными ферменного типа (рис. 1 в). 


Основной 
элемент 


[В 
Встроенный 
элемент 
а) 6) 


в) 


Рис. 1. Схемы моделирования армирования: а — дискретно-распределенная; 
б — со встроенными стержневыми элементами; в — основным и встроенным элементами (рисунок авторов) 


В российской практике прочностного расчета монолитных зданий и сооружений из железобетона преиму- 
щественно используются узко специализированные программные комплексы ЛИРА-САПР и $САП- 
ОЁЯсе [4, 5], в которых для моделирования рамно-связевых каркасов применяется технология ансамблирования 
пластинчатых и балочных КЭ. При этом в зонах сопряжения перекрытий и колонн автоматически вводятся так 
называемые жесткие вставки, представляющие собой звездообразно расположенные балочные КЭ с искус- 
ственно завышенной изгибной жесткостью (рис. 2). Геометрия жесткой вставки соответствует поперечным 
размерам сечения колонны. Такой подход позволяет придать более физичный характер распределению изгиба- 
ющих моментов в КЭ перекрытия. 


Балочные КЭ 
Пластинчатые КЭ 


* ие ар. 
аи 
ых 
Балочные КЭ с модулем 
упругости 10*.Е 


а) 6) 


Рис. 2. Моделирование зоны сопряжения перекрытия и колонны: 
а — зона сопряжения; 6 — расчетная схема зоны сопряжения с жесткой вставкой (рисунок авторов) 


Альтернативой жестким вставкам с балочными КЭ с завышенной жесткостью является процедура коопти- 
рования (связывания) степеней свободы узлов пластинчатых элементов, примыкающих к узлу стержневого 
элемента, с соответствующими узловыми перемещениями и углами поворотов стержневых КЭ. Естественно, 
при таком подходе необходимо в местах сопряжений пластинчатых и стержневых КЭ предусмотреть сгущение 
сетки до размеров поперечных сечений колонн. Отметим, что рассмотренный подход также может быть про- 
граммно автоматизирован. 

Материалы и методы. В настоящее время разработчики конечно-элементных программных комплексов 
при построении матриц жесткости пластинчатых и оболочечных конечных элементов широко используют 
алгоритм МГТС (М хе4 Пиегро]аНоп оЁ Тепзопа1 Сотропеп($), базирующийся на процедуре независимой (раз- 
дельной) аппроксимации изгибных и сдвиговых деформаций. Цель данной процедуры состоит в исключении 
эффекта «заклинивания» или ложного сдвига. 

Для моделирования колонн монолитных железобетонных каркасов применяют прямолинейные двухузловые 


балочные КЭ с шестью степенями свободы в узле, которые включают три перемещения в направлении локаль- 
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ных осей и соответствующие угловые перемещения. 

В программном комплексе АМЗУ$ Месвашса| для расчета пространственных рамно-связевых каркасов 
предусмотрены пластинчатый четырехузловой КЭ $НЕЧ.63 и двухузловой балочный КЭ ВЕАМ!88'. Анализ 
трехмерного напряженно-деформированного состояния конструкций из железобетона выполняют с помощью 
специального восьмиузлового элемента ЗОГ.О65. 

Результаты исследования. Выполним расчет напряженно-деформированного состояния монолитного трех- 
этажного железобетонного каркаса с шагом колонн 5х7 м. Высота этажа — 4,7 м; размер сечения колонны — 
0,4х0,4 м; толщина перекрытия — 0,2 м. Модуль упругости, коэффициент Пуассона, удельный вес материалов 
перекрытий Е = 27510“ МПа; 1 = 0,2; у = 2440 ким 3, колонн Ё = 3х10% МПа; и = 0,2; у = 2500 ким 3. Расчетные 


значения сопротивления бетона осевому сжатию и растяжению соответственно составляют: В, = 25,5 МПа, 
В, =2,37 МПа. 


Полагаем, что на все перекрытия каркаса действует постоянная равномерно распределенная нагрузка 
Ч= 2 кПа. Расчеты выполняем с учетом собственного веса каркаса. Считаем, что основания колонн первого 
этажа каркаса жестко закреплены. 

На первом этапе расчета перекрытия и колонны каркаса моделируем соответственно пластинчатыми 
ЗНЕГ.63 и балочными ВЕАМ188 КЭ. 

В таблице 1 приведены результаты сравнительных расчетов перекрытия первого этажа каркаса. Вычисления 
проводились с шагом сетки на перекрытиях 0,2 м и 0,5 м. Колонны в обоих случаях разбивались на 6 КЭ. В 


таблице 1и далее обозначено: и, — максимальный прогиб; М», м === изгибающие моменты относительно 


глобальных осей. В третьей строке таблицы 1 приведены данные расчета каркаса с жесткими вставками (шаг 
сетки 0,2 м). Жесткие вставки с размерами сторон 0,4х0,4 м моделировались пластинчатыми КЭ ЗНЕ 1.63 с 
модулем упругости 10 3х. 

Отметим, что значения М, тах, Му тах, Приведенные в таблице 1, относятся к небольшим локальным зонам 
сопряжения колонн и перекрытий. Величины М; ть и М, ты распространяются по периметру каркаса между 
рядами крайних колонн. 

Аналогичный расчет каркаса был выполнен с использованием программного комплекса ЛИРА-САПР. Зна- 
чение максимального прогиба на первом этаже при аналогичном нагружении и механических постоянных ма- 


териала, полученное с помощью программного комплекса ЛИРА САПР, составило и, = -3,06 мм, что сопоста- 
* 
вимо с расчетом в АМЗУ$ и, = -3,49 мм при шаге сетки на перекрытиях 0,2 (таблица 1). При расчете с 


использованием комплекса ЛИРА-САПР на шаге построения «аналитической модели» в точках пересечения 
колонн и перекрытий включались «контуры продавливания», опосредованно влияющие на изгибную жесткость 
перекрытий в сторону ее увеличения. 
Таблица 1 
Значения экстремальных изгибающих моментов в зоне сопряжения перекрытия первого этажа и колонны 
при использовании НЕТ Т.63 


Шаг сетки, Число вы Му, кНм Му, кНм 
м неизвестных тт ри тт р 
0,5 103032 —0,004908 —13,1 37,2 —20,6 45,3 
0,2 583200 —0,005409 —15,2 70,2 —22,5 99,7 
0,2" 583200 —0,003492 —10,8 106 —17,5 106 


В таблице 2 приведены значения экстремальных изгибающих моментов для перекрытия первого этажа, 

* 
полученные с помощью АМ$ЗУ$ (шаг сетки 0,2 м) и ЛИРА-САПР (шаг сетки 0,395 м). Следует отметить, что в 
комплексе ЛИРА-САПР для количественной оценки значений М, и м введена шкала погонных изгибающих 
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моментов, т. е. результирующие моменты приведены к полосе плиты шириной 1м. 


' Басов К. А. АМ$У5: справочник пользователя. М.: ДМК Пресс, 2005. 640 с. 309 
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Таблица 2 
Сравнение экстремальных значений изгибающих моментов в зоне сопряжения перекрытия первого этажа 
и колонны, вычисленных с помощью АМЗУ$ и ЛИРА-САПР 


‚ Му, кН.м Му, кНм 
Программный комплекс - - 
тт тах тт тах 
АМ$У$ —10,8 106 =, 5 106 
ЛИРА-САПР —31,5 8,37 —42,2 12,9 


Из представленных в таблице 2 данных видно, что полученные с использованием АМЗУ$ и ЛИРА-САПР 
одноименные экстремальные моменты существенно различаются по величине. Данное обстоятельство объясня- 
ется наличием в исследуемой конечно-элементной модели локальных зон концентрации внутренних усилий в 
местах сопряжения пластинчатых и стержневых КЭ. Причем принятые размеры пластинчатых КЭ в обоих ком- 
плексах не позволяют точно смоделировать градиенты изменения изгибающих моментов в указанных зонах 
концентрации. 

Анализируя результаты расчета каркаса по пластинчато-стержневой схеме приходим к выводу, что данная 
модель не позволяет исследовать детальную картину напряженного состояния конструкции. В частности, не 
удается проанализировать зоны растягивающих нормальных напряжений 0х, 0, 0», возникающих в местах со- 
пряжения колонн и перекрытий. В этой связи актуальной является задача построения расчетной модели повто- 
ряющегося фрагмента каркаса, позволяющей выполнить численное исследование объемного напряженно- 
деформированного состояния, в том числе с учетом армирования. На рис. 3 изображен план перекрытия рас- 
сматриваемого каркаса с выделенными повторяющимися фрагментами а, Ь, с. Область повторяющегося фраг- 
мента а является объектом дальнейшего исследования. 


Рис. 3. Повторяющиеся фрагменты а, Ь, с перекрытия каркаса (рисунок авторов) 


Для моделирования объемного напряженно-деформированного состояния повторяющегося фрагмента при- 
меним восьмиузловой КЭ типа $ОГ.ТО185 с тремя степенями свободы в узле. 

На рис. 4 представлена расчетная схема повторяющегося фрагмента а рассматриваемого каркаса. Приведен- 
ная расчетная схема относится к первому этажу каркаса. Здесь сосредоточенная сила Р =72,6 кН взята из эпюры 
продольных сил, полученной с использованием пластинчато-стержневой модели каркаса. При этом сосредото- 
ченную силу Р преобразуем к статически эквивалентному давлению 4, = 453,5 кПа, действующему на площад- 


ке 0,2х0,2 м (14 часть сечения колонны). Статические граничные условия Чх, Чу, И, накладываются на узлы 


КЭ с учетом циклической симметрии картины деформации перекрытия. Соответствующая конечно-элементная 
модель фрагмента, построенная на базе объемных КЭ типа $ОГ./1О185, приведена на рис. 5. В данном случае 
шаг сетки объемных КЭ принят равным 0,1 м. Давление на фрагмент перекрытия и 1/4 части колонны прикла- 
дываем с помощью элемента $0ВЕ1 542. 

Результаты сравнения жесткостных свойств пластинчато-стержневой и объемной конечно-элементных мо- 


делей рассматриваемого фрагмента каркаса в виде картин распределения перемещений Ч, приведены на 


рис. 6. Картины распределения Чу на рис. 6 б и рис. 6 в соответствуют расчетам фрагмента с учетом и без уче- 


та силы Р. Как видно из представленных результатов, расчет без учета силы Р (рис. 6 в) дает значения переме- 


2 То же. 
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щений, наиболее близкие к данным пластинчато-стержневой модели. Дальнейшие численные эксперименты 
будем проводить с учетом нагружения колонны фрагмента давлением 4,. 


и. =0 ры и 
БИ - их =0 
у и. 
Х 
и, =0 
их 


и, =0 


и, =0 
и, =0 
Е Рис. 5. Конечно-элементная модель фрагмента 
Рис. 4. Расчетная схема фрагмента (рисунок авторов) (рисунок авторов) 
-4 : -3 н -3 
-0,757.10 ^ 0,761.10 2 - -0,228-10 


-0,3603-10` -0,3098-10`2 


-0,3112.10-2 


-2 
5 -0,3757-10 
-0,3279-10 -0.1984-1072 -0,1478-1072 


а) 6) в) 


«: -0,3252-1072 
-0,1471-10 


Рис. 6. Результаты вычисления перемещений и, по двум моделям: а — пластинчато-стержневая модель; 


б — объемная модель с учетом силы Р; в — объемная модель без учета силы Р (рисунок авторов) 


Выполним исследование влияния армирования на напряженно-деформированное состояние рассматривае- 
мого фрагмента каркаса. Считаем, что перекрытия каркаса безригельные и бескапительные. Рассмотрим три 
схемы моделирования армирования части перекрытия, принадлежащей фрагменту: | — дискретное армирова- 
ние с использованием стержневых КЭ (ВЕАМ188); 2 — дискретное армирование с использованием ферменных 
КЭ (ЛМК 180); 3 — распределенное армирование с использованием пластинчатых КЭ (ЗНЕГЛ 63). Следует от- 
метить, что моделирование арматуры посредством ферменных КЭ также подходит для расчета преднапряжен- 
ных железобетонных конструкций? [1, 6, 7]. 

Конструктивные элементы дискретного армирования фрагмента каркаса продемонстрированы на рис. 7. 
Материал армирующих стержней — сталь (Е = 2х10° МПа; и = 0,28; у = 7800 км °). 

Армирующие сетки выполнены из стержней следующего диаметра: фоновое армирование (по всей плоско- 
сти перекрытия) — 10 мм; армирование ригелей — 12 мм; армирование капители — 16 мм. Диаметры стержней 
для поперечного армирования ригелей и капителей имеют значения, аналогичные вышеприведенным. Армиру- 


ющие стержни колонны диаметром 12 мм расположены по углам сечения. 


300 


3500 3500 


200 
2500 


а) 6) ы 


Рис. 7. Конструктивные элементы дискретного армирования: а — фоновое верхнее/ нижнее армирование; 


б — армирование ригелей; в — армирование капители (рисунок авторов) 


3СП 52-101-2003. Бетонные и железобетонные конструкции без предварительного напряжения арматуры. М., 2004. 59 с. 
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Назначенные значения диаметров армирующих стержней взяты на основании прототипа [8—10]. Для моде- 
лирования арматуры используем балочные или ферменные КЭ. 

Альтернативой модели дискретного армирования является подход, базирующийся на принципе «размазыва- 
ния», так называемая модель распределенного армирования. Суть принципа «размазывания» в общем случае 
состоит в введении приведенного модуля упругости, являющегося функцией от модулей упругости компонен- 
тов неоднородного материала и их объемных концентраций. Однако расчетная практика [1] показала, что для 
железобетонных конструкций, работающих на изгиб, более реалистичные результаты дает введение дискрет- 
ных армирующих слоев, эквивалентных объему, занимаемому армирующими стержнями. Геометрия таких сло- 
ев повторяет геометрию армирующих сеток. Для рассматриваемого фрагмента каркаса на рис. 8 показаны экви- 
валентные армирующие слои, моделируемые пластинчатыми КЭ типа ЗНЕТ 1.63. 


6) в) 
Рис. 8. Эквивалентные армирующие слои: а — слой для фонового армирования (верхнего/ нижнего); 
б — слои армирования ригелей; в — слои армирования капители (рисунок авторов) 


На основании условия равенства объемов армирующих элементов и совпадения геометрии зон армирования 
получим следующие толщины для эквивалентных армирующих слоев: слой фонового армирования — 1,6 мм; 
слои армирования ригелей — 1,1 мм, 0,74 мм, 1,2 мм, 0,78 мм; слои армирования капители — 5 мм, 2,7 мм. 


В процессе вычислительного эксперимента исследовалась зависимость прогиба и, в точке К и максималь- 


й о, и: от принятой схемы армирования исследуемого фрагмен- 


у тах ? 2 тах 


ных растягивающих напряжений о 


х тах ?› 


[@} 


та (рис. 4). В таблице 3 приведены результаты расчетов. 


Таблица 3 
Значения прогибов и напряжений при различных способах моделирования арматуры 
ИОДЕЛирОВаНИЕ и, в точке К, мм о „› МПа 0’ м„› МПа о, их> МПа 
арматуры 
Без арматуры —3,32 3,26 4,40 1,44 
ВЕАМ188 —2,81 2,77 3,76 1,08 
ГМК 180 —2,92 2,75 3,74 1,06 
$НЕ!Г 63 —2,06 1,19 1,74 0,513 


Анализ результатов показал, что применение пластинчатых КЭ приводит к существенному занижению мак- 
симальных растягивающих напряжений. Данные, полученные с использованием балочных и ферменных КЭ, 
практически совпадают, что объясняется низкой изгибной жесткостью армирующих стержней. 


Важно подчеркнуть, что для рассмотренного варианта нагружения при использовании стержневых КЭ 


(ВЕАМ188 и 11МК180) в зоне сопряжения колонны и перекрытия условие прочности для напряжений ок и 


х тах 


Иа 


[@} не выполняется. 


у тах 
На рис. 9 приведена визуализация процесса трещинообразования в плите перекрытия, полученная при моде- 
лировании бетона объемными восьмиузловыми КЭ $О0110565. Следует отметить, что в перекрытии появляются 


только микротрещины, граница которых хорошо согласуется с полем напряжений сдвига [ер в плоскости пе- 


рекрытия (рис. 10). 
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Рис. 9. Визуализация процесса трещинообразования в плите перекрытия: а — общий вид фрагмента; 
б — вид зоны трещинообразования сверху; в — вид зоны трещинообразования снизу (рисунок авторов) 
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Рис. 10. Картина распределения напряжений о’„, в перекрытии: 


а — вид фрагмента сверху; 6 — вид фрагмента снизу; в — шкала напряжений (рисунок авторов) 


Для моделирования физической нелинейности бетона использовалась модель Уиллама-Варнке [11] со сле- 
дующими параметрами: коэффициент передачи сдвиговых усилий при открытой трещине — 0,3; коэффициент 
передачи сдвиговых усилий при закрытой трещине — 0,7; коэффициент понижения жесткости при образовании 
трещины в результате растяжения — 0,6. 

Обсуждение и заключения 

1. На основании выполненного линейно-упругого расчета монолитного каркаса из железобетона доказано, 
что пластинчато-стержневая модель, широко используемая в практике проектирования, не позволяет количе- 
ственно оценить величины растягивающих нормальных напряжений в зонах сопряжения плит перекрытий и 
КОЛОНН. 

2. Разработана методика численного моделирования объемного напряженно-деформированного состояния в 
месте сопряжения перекрытия и колонны монолитного железобетонного каркаса, позволяющая оценить факти- 
ческий запас прочности данного узла, а также уточнить несущую способность соответствующего здания или 
сооружения при различных сценариях нагружения. 
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